Gravitational collapse of an isentropic perfect fluid with a linear equation of state 

Rituparno Goswamjj and Pankaj S JoshjJ 
Department of Astronomy and Astrophysics 

Tata Institute of Fundamental Research 
Homi Bhabha Road, Mumbai 400 005, India 

We investigate here the gravitational collapse end states for a spherically symmetric perfect fluid 
with an equation of state p = kp. It is shown that given a regular initial data in terms of the 
density and pressure profiles at the initial epoch from which the collapse develops, the black hole or 
naked singularity outcomes depend on the choice of rest of the free functions available, such as the 
velocities of the collapsing shells, and the dynamical evolutions as allowed by Einstein equations. 
This clarifies the role that equation of state and initial data play towards determining the final fate 
of gravitational collapse. 
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Considerable interest is seen in recent years to examine the final fate of gravitational collapse of a massive matter 

cloud within the framework of Einstein's theory of gravity. This is due to the importance of this problem in black 

hole physics, and the related cosmic censorship conjecture which states the spacetime singularities of collapse must 

^ . be hidden within black holes, and would not communicate with faraway observers in the spacetime (for some recent 

^ ■ reviews, see e.g. [1-5]. The assumption that continual collapse of a matter cloud, such as a massive star which 

^^ , has exhausted its nuclear fuel, ends only in a black hole (BH) and not naked singularity (NS) is crucial to many 

of the considerations in the physics and astrophysics of black holes. Several dynamical collapse scenarios have been 

investigated extensively from such a perspective, which include the radiation collapse within the framework of a Vaidya 

(^ ■ metric (see e.g. ,6], and references therein), collapse of a dust cloud [6-10], some considerations on perfect fluids (many 

T^lj- ' of these being numerical) [11-19], massless scalar fields [20-26], and also more general matter fields [27-29]. 

C^ ■ Our purpose in this note is to study analytic models of spherical gravitational collapse of a perfect fiuid with an 

^j ■ equation of state p = kp. This is of interest because this is a well-studied and extensively used case in astrophysics, 

^^i which offers a physically interesting model. While modeling a realistic collapse, one has to consider the equation of 

' ' state (EOS) of the collapsing matter as an additional constraint over the Einstein equations. This is a long standing 

j^i interesting question in the arena of dynamical collapse theories as to how a physically realistic EOS affects the 

!* ■ evolution of a collapse in terms of it's final state. In the later stages of collapse one could not ignore the pressures, 

. ^ ' and hence the case when matter is an isentropic perfect fluid with a linear equation of state offers a useful scenario 

S^ r to examine collapse in ultra- relativistic limits. 

5h ' We show here that a perfect fluid collapse could end in cither of the BH/NS flnal states, depending on the nature 

of the initial data, and the allowed evolutions for the collapsing matter. Given a regular initial data for matter in 
terms of the regular density and pressure proflles, the exact outcome in terms of the above depends on the choice of 
rest of the free functions available, such as the velocities of the collapsing shells, and the allowed evolutions for the 
collapsing matter. Our results thus provide some insight on the role that equation of state and initial data play to 
determine the final fate of continual gravitational collapse in terms of the BH/NS end states. 

The basic regularity conditions and the Einstein equations are given in Section 2. The collapse is examined in 
Section 3, and in Section 4 we study how the nature of the singularity is determined by the initial data and the 
evolution chosen. Some conclusions are summarized in Section 5. 
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II. EINSTEIN EQUATIONS AND REGULARITY CONDITIONS 

The spacetime geometry within the spherically synimetric collapsing cloud can be described by the metric in the 
comoving coordinates (i, r, 9, 0) as given by, 

ds2 = -e^-i'-^Ut' + e2^(*''-)dr2 + R\t, r)dn^ (1) 

where dft^ is the line element on a two-sphere. In this reference frame the energy-momentum tensor for a perfect 
fluid is given by, 

T'^ = {p+p)V'V^ +pg'^ (2) 

where V^ is unit timelike vector. We also take the matter fields to satisfy the weak energy condition^ i.e. the energy 
density measured by any local observer is non-negative. Then we must have, 

T.kVW' > (3) 

which amounts to, 

P>0;p + P>0; (4) 

Since we consider here an isentropic fluid, whose pressure is a linear function of the density only, the equation of state 
of the collapsing matter is given by, 

p{t,r) ^ kp{t,r) (5) 

where k e (—1, 1] is a constant. We note that the cases of negative k are the dark energy fluids with negative pressures. 
Now for the metric ^ the Einstein equations take the form (in the units SttG — c = 1) 

F' 1 F 

(6) 



R'^R' k mk 

k 



V = — 



/c+1 



[ln(p)]' (7) 



R'G - 2Rv'G = (8) 

G-i/ = l-| (9) 

where, 

G(i, r) = e-2'>(i?')2; H{t, r) ^ e-^" {Rf (10) 

The arbitrary function F = F{t, r) here has an interpretation of the mass function for the cloud, and it gives the total 
mass in a shell of comoving radius r on any spacelike slice t — const. We have _F > from the energy conditions. In 
order to preserve the regularity at the initial epoch, we have F(ti,0) = 0, that is, the mass function should vanish 
at the center of the cloud. Since we are considering collapse, we have R < 0, i.e. the physical radius R of the cloud 
keeps decreasing in time and ultimately reaches R — 0, which denotes the singularity where all matter shells collapse 
to a zero physical radius. We can use the scaling freedom available for the radial co-ordinate r to write R = r a.t the 
initial epoch t = U. Let us introduce a function v{t, r) as defined by, 

v{t,r) = R/r (11) 

we then have i?(i, r) = rv{t, r), and 

v{U,r)^l; v{ts{r),r) =0; v<0 (12) 



The time t — ts{r) corresponds to the shell- focusing singularity at i? = 0, where all the matter shells collapse to a 
vanishing physical radius. 

The description of the singularity in terms of v{t, r) has the following advantage. The physical radius goes to the 
value zero at the shell- focusing singularity, but we also have i? = at the regular center of the cloud at r = 0. 
This is to be distinguished from the genuine singularity by the fact, for example, that the density and other physical 
quantities including the curvature scalars are all finite at a regular center r = 0, even though i? = holds there. This 
is achieved, as we point out below, by a suitable behaviour of the mass function, which should go to a vanishing value 
sufhciently fast in the limit of approach to the regular center where even though R goes to zero the density must 
remain finite. On the other hand, we note that at t = U we have v = 1, and then as the collapse evolves v continuously 
decreases to become zero only at the singularity, i.e. v = uniquely corresponds to the genuine spacetime singularity 
at i? = 0. 

We thus see that there are now five total field equations with five unknowns, which are p, ip, v, R, and F. Solutions 
of these equations, subject to the weak energy condition and the given regular initial data for collapse at the initial 
spacelike surface t = ti, determine the matter distribution and metric of the space-time. We then have specific time 
evolutions of the initial data which define the collapse final states. It turns out that there exist classes of solutions, 
which give either a black hole or a naked singularity as the end state of the collapse, depending on the nature of the 
initial data, and the class of evolutions chosen. 

III. COLLAPSING MATTER CLOUDS 

It is now possible to consider the gravitational collapse of a perfect fluid within this framework as we discuss below. 
The regularity conditions as defined above set up the initial data at the initial surface t = ti from which the collapse 
develops. We now consider the collapse equations which allow us to see when the singularity occurs, and how the 
initial data and the classes of evolutions as governed by the Einstein equations lead to the formation of the spacetime 
singularity. 

We consider a general mass function F{t,r) for the collapsing cloud, which can be written as 

F{t,r) =r^Mir,v) (13) 

where A^ is a regular and suitably differentiable function, and A^ > 0. As seen from equations below, the regularity 
and finiteness of the density profile at the initial epoch t ^ ti, and at all other regular epochs before the singularity at 
R = develops, requires that F goes as r^ close to the center. Hence we note that since A^ is a general (at least C^) 
function, the equation equation (|13|l is not really any ansatz, or a specific choice, but quite a generic form of the mass 
profile for the collapsing cloud. We note that [19] considered mass profiles which are analytic (with F{t, r) — F{r, R)), 
and perfect fluids could form naked singularity, but we do not impose such an assumption on mass function here. 
Then equation IJBJ gives, 

3M+r[Mr+M^.v'] 1 M „ 

P= 2?—, 7T~ = ^7^2- (14) 

v^yu + rv ) k v 

The regular density distribution at the initial epoch is given by, 

Po{r)^iM{r,l)+rM{r,l)^r (15) 

It is evident that, in general, as w — > 0, p — > cx). That is, both the density and pressure blow up at the singularity. 
We note from equation 114() that p = p{r, v) and hence v' = f{r, v). Now rewriting equation (|14() we get, 

SkM + krM.r+Q{r,v)M^y^O (16) 

where, 

Q{r,v) = (k + l)rv' + v (17) 

Now the above equation (|16|) has a general solution of the form [S^l , 

T{X,Y) = (18) 

where X{r,v,A4) and Y{r,v,A4) are the solutions of the system of equations, 

dAi dr dv 
3k kr Q 



Amongst all the classes of solutions of ^A (r, v) as given by equation (|18|l , only those are to be considered which obey 
the energy condition, which are regular, and which in the limit of w — > give p —^ oo. In other words, the equation of 
state as given by the perfect fluid condition p = kp, and the energy condition isolate the class of the mass functions 
to be considered. 

We can now directly integrate equation lO to get, 

Hr,v) = -^MP)] (20) 

Let us now define a suitably differentiable function A{r, v) in the following way, 

iy\r,v)^Air,v)^,R' (21) 

That is A(r,v)^v = v' / R' . We note that equation (14) can in principle give solutions which are not regular, or such 
that the function A{r, v) as defined above is not regular. However, regularity of A4 and j4 is a necessary assumption 
we make here, which will be used later also (see e.g. remarks after equation (30)). Hence, only the class of regular 
solutions is to be considered. 

Our main interest here is in studying the shell-focusing singularity at i? = which is physical singularity where 
all the shells collapse to zero radius. Hence we assume that there are no shell-crossing singularities in the spacetime, 
where R' = and so the function A{r, v) is well-defined. 

Some comments are in order here on our assumption that R' > 0, that is, we consider the situation with no shell- 
crossing singularities. This is because, it is generally believed (see e.g. l3l| that such singularities can be possibly 
removed from the spacetime as they are typically gravitationally weak, and also because spacetime extensions have 
been constructed through the same in certain cases. Under the situation, we are interested only in examining the 
nature of the shell-focusing singularities at i? = 0, which are genuine curvature singularities, where the physical radii 
for all collapsing shells vanish, and the spacetime necessarily terminates without extension. 

Specifically, R' > implies that we must have v + rv' > 0. Since v is necessarily positive, it follow that this will 
be satisfied whenever v' is greater or equal to zero, or even when it is negative the magnitude of rv' should be less 
then that of v. Later in this section we shall derive an expression for the quantity v' , in terms of the initial data 
and the other free evolutions as allowed by the Einstein equations. Hence it follows that we can specifically state the 
condition for avoidance of shell-crossings in terms of the behaviour of these functions. In particular, it turns out that 
whenever the singularity curve tg (r) (which corresponds to R ~ 0) is increasing (or when it decreases at a sufficiently 
slow rate) with a slope greater or equal to zero at the origin, the shell-crossing singularities are avoided at least in 
the vicinity of the regular center r — 0. We then have a ball of finite radius around the regular center which contains 
no shell-crossings till the final singularity formation at _R = 0. We shall, however, not go into further details here. 

At the initial epoch we have. 



Mr,v)^v\,^^ = --r— 



P'oJr) 
1 [po{r) 



(22) 






In fact, for all epochs, the relation between the function A4 and A is given by equation 10 as A^^R' = — -j^ In 
If we consider a smooth initial profile, i.e the gradient of the initial density vanishes at the center, then we must have 
A{r,v) = rg{r,v), where g{r,v) is another suitably differentiable function. 
Now using equation (|21|) we can integrate (jS} to get, 

G(r, w) = 6(r)e2'^^ (23) 

Here h{r) is another arbitrary function of the comoving coordinate r. A comparison with dust collapse models 
interprets b{r) as the velocity function for the collapsing shells. Following this parallel, we can write, 

h{r) = \ + r%o{r) (24) 

Finally, using equations (f^ , (EOJl and if^ in (0 we have. 



\fvv = -p ''+^Je'^''^vbQ{r)+vh{r,v)+M{r,v) (25) 

where, 

h{r,v) = ^^^ (26) 



Integrating the above equation we have, 

t{v,r)^ —, ^ (27) 

Jv p ''+1 ^/e^^'^vbo + vh + M 

Note that the variable r is treated as a constant in the above equation. Close to the center we can write t{v, r) as, 

t{v,r)^t{v,0) + rX{v)+O{r^) (28) 



Here the function X{v) is given by. 



where, 



Xi.) = -U\v^^^im (29) 



B{r, v)^ p~T^ y/e'^'-^vbo+vh + M; Bi ^ B^r (30) 

In order to obtain equation (28), we note that we require the integral (27) could be differentiated. This is possible 
because it is finite by definition, and then we need all the functions, namely A(r,v),bo{r) and Al{r,v) to be suitably 
differentiable. In our case we require them to be at least C^ for r not equal to zero, and C^ for r = 0. 

Thus we see that the time taken for the central shell to reach the singularity is given as 

Jo B{0,v) 

The time for other shells to reach the singularity is given by the following, which defines the singularity curve 
developing in the spacetime as end result of collapse, 

i,(r)-i,„+rA'(0) + O(r2) (32) 

It is now clear that the value of the quantity A'(O), which represents the tangent to the singularity curve, depends 
on the functions 6o, -M and h, which have the initial values as dictated by the initial data at t — ti, and which are 
functions of r and v as the case may be. Hence, a given set of density and velocity profiles, together with the evolutions 
chosen, completely determines the tangent at the center to the singularity curve. Further, from equation (|25|l . we get, 
for a constant v surface, 

V^v ^X{v)B{0,v)+O{r) (33) 

We note that as seen above, X{0) involves functions depending on the initial data, and also the evolutions A and 
A4. The relation between the functions A and M is given by the perfect fiuid equation of state, for which we have 
shown that solutions exist, and the Einstein equations as noted earlier. Among different classes of solutions only those 
are to be considered which ensure the density and pressures to blow up at the singularity. 

One has now to understand the structure of this singularity, and to examine when it will be visible, and when 
covered within an event horizon of gravity, i.e. hidden within a black hole. 

IV. NATURE OF THE SINGULARITY 

It is now possible to see for a perfect fluid collapse with a linear equation of state, how the initial data and the 
allowed evolutions determine the final fate of collapse in terms of either a black hole or a naked singularity. The 
apparent horizon within the collapsing cloud is given hy R = F. If the neighborhood of the center gets trapped 
earlier than the singularity, then it is covered and a black hole results, otherwise it is visible with non-spacelike future 
directed trajectories escaping from it. In other words, we examine below when there will be families of null geodesies 
existing, which will be future directed and outgoing, and which terminate in the past at the singularity, thus making 
the communication from the singularity to an outside observer possible, as opposed to a black hole situation where 
this will not be the case. 

In order to consider the possible existence of such trajectories and to examine the nature of the central singularity 
at i? — 0, r = 0, let us consider the equation for outgoing radial null geodesies which is given by, 

^ = e'^-' (34) 

dr 



The singularity occurs at v{ts{r),r) = 0, i.e. R{ts{r),r) — 0. Therefore, if there are any future directed null geodesies 
existing, which terminate in the past at the singularity, we must have i? — > as i — > ts along these curves. Now 
writing equation H^M() in terms of variables (u = r"', R) , we have. 



du a 



'-§''-' 



Choosing a = | and using equation @ we get. 



(35) 



(36) 



If there are null geodesies which terminate at the singularity in the past with a definite tangent, then at the 
singularity we have ^ > 0, in the {u, R) plane with a finite value. Hence it follows that all points r > on the 




singularity curve are covered necessarily, because F/ R — > oo with ^ ^ — cx) for any of these, and hence no outgoing 
null geodesies can terminate at these points in the past. The central singularity at r = could however be naked. 
Define the tangent to the outgoing null geodesic from the singularity as. 



R 



Xn = lim lim — 



dR 
du 



Using equation (jSS)l and (|S5|) . we then get. 



xi - -/b(o:o)a'(o) 



(37) 



(38) 



Let us now deduce the necessary and sufficient conditions for a naked singularity to exist, that is, for null geodesies 
with a well-defined tangent to come out from the central singularity. Suppose we have X{Q) > 0, then we always have 
(from equation (|38(l ). a;o > and then in the (_R, u) plane, the equation for the null geodesic that comes out from the 
singularity is given by 



R 



XqU 



(39) 



In other words, equation 139|) is a solution of the null geodesic equation in the limit of the central singularity. Thus 
given X(0) > 0, we can always construct a solution of radially outgoing null geodesies emerging from the singularity. 
This makes the central singularity visible. In the (t, r) plane, the null geodesies outgoing from the singularity will be 
given as. 



^-^^(O) = xor3 



(40) 



It follows that X{0) > implies xq > and we get radially outgoing null geodesies emerging from the singularity, 
giving rise to the central naked singularity. 

On the other hand, if X{Q) < 0, then we see that the singularity curve is a decreasing function of r. Hence the 
region around the center gets singular before the central shell, and the spacetime then terminates there. In this case, 
if there were any outgoing null geodesic from the central singularity, it must then go to a singular region, or outside 
the spacetime which is impossible. Hence when X{0) < 0, we always have a black hole solution. 

If A'(O) = then we will have to take into account the next higher order non-zero term in the singularity curve 
equation, and do a similar analysis by choosing a different value of a in equation (|35|) . 

We have thus shown above that X{0) > is the necessary and sufficient condition for null geodesies to come out 
from the central singularity with a definite positive tangent. It is thus seen how the initial data, together with the 
evolutions chosen in terms of the free functions such as 6o, A^ and h, fully determine the final end product of collapse 
in terms of either a black hole or a naked singularity. This is as determined by the X{0) values above, or more generally 
in terms of the behaviour of the singularity curve in the vicinity of the central singularity. This is because X(0) is 
determined by these initial profiles and the evolutions chosen as given by equation H29|l . which in turn determine the 
end states. Therefore, given any initial regular density and pressure profiles for the matter cloud from which the 
collapse develops, there always exist velocity profiles for collapsing matter shells, and evolutions as determined by the 
Einstein equations, so that the end state of the collapse would be BH or NS, depending on the choice made. 

As seen above, the different outcomes are characterized by the positive or negative values of X(0). This in turn 
depends on the initial values of the functions such as the initial mass and velocity profiles, and the evolutions chosen, 
as permitted by the Einstein equations. Thus the measures of outcomes leading to BH/NS phases are accordingly 
decided (e.g. the entire set of initial data and evolutions giving X{0) > leads to NS). 



V. CONCLUSIONS 

While some of the numerical models for perfect fluids indicated that naked singularities could arise for only a 'soft' 
equation of state, our results point out that within a generic perfect fluid collapse scenario with equation of state 
p = kp, as such the value of k chosen does not appear to have any special significance. What matters is the initial 
data, and the chosen evolutions (i.e. the classes of allowed solutions to the Einstein equations), which then take 
this given initial data to a specific outcome, depending on the choice made. Also, if for a given chosen evolution, if 
the value X{0) was negative (or positive), then such will be the case by continuity for all neighbouring or close by 
evolutions where 'nearness' is defined in some suitable sense. Hence these outcomes in terms of BH or NS may be 
considered stable in a certain sense (within spherical symmetry) as characterized above. We thus see that there are 
classes of solutions to Einstein equations for perfect fluid models where given the matter initial data at the initial 
surface t — U, these evolutions take the collapse to end up either as a black hole or the naked singularity, depending 
on the choice of the class. This also means that the total space of evolutions can be divided into distinct subspaces, 
those that evolve a given initial data into black holes, and others that go to a naked singularity. The results on dust 
collapse are of course contained here as a special case with p — 0. 

We should mention here that though the above give some information on dynamical evolution of collapse while 
pressures are included in the analysis, non-spherical perturbations will have important inputs to decide on the issues 
such as genericity and stability. It is possible that methods such as those developed in [32-33] could be useful in that 
direction. 
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